Decoupling and Depinning II: Flux lattices in disordered layered superconductors 
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Phase transitions of a flux lattice in layered superconductors with magnetic field perpendicular 
' to the layers and in presence of disorder are studied. We find that disorder generates a random 

Josephson coupling between layers which leads to a Josephson glass (JG) phase at low temperatures; 
vanishing of the JG order identifies a depinning transition. We also find that disorder and thermal 
fluctuations lead to layer decoupling where the renormalized Josephson coupling vanishes. Near 
decoupling an anharmonic regime is found, where usual elasticity and the resulting Bragg glass are 
, not valid. The depinning line crosses the decoupling line at a multicritical point, resulting in four 

I ■ transition lines. The phase diagram is consistent with the unusual data on BiiSr-iCaCuiO% such 

' as the "second peak" and depinning transitions. The Josephson plasma frequency is evaluated in 

the various phases. 

PACS numbers: 74.25.Qt,74.25.Dw,74,50.+r 
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■ I. INTRODUCTION 

g : 

c/2 ' The phase diagram of layered superconductors in a magnetic field B perpendicular to the layers is of considerable in- 
terest in view of extensive experiments on high temperature superconductors*^. A first order transition in Y Ba^Cu-^O'j 
(YBCO) and in Bi2Sr2CaCuiO% (BSCCO) has been interpreted as a melting transition of the flux lattice. The data 
suggests that the first order line terminates at a multicritical point, which for BSCCC^ 3 - is at Bq ss 300 — 10 3 G and 
T « 40 - 50K, while for YBCQi it is at Bq « 2 - 10T and Tq w 60 - 80-ftT, depending on disorder and oxygen 
concentration. The multicritical point also terminates a "second peak" transitions* 2 *^ which is manifested by a sharp 
increase in magnetization and in critical current. The transition line at B ss Bq and T < Tq is weakly T dependent 
and was found, for BSCCO, to be smoothly connected with the first order line-V Neutron scattering and /iSR datai*& 
show that positional correlations of the flux lattice are significantly reduced near these phase boundaries, except 
' however, near the multicritical point where a reentrant behavior is observed^. Data on Ndi^Ceo.i^CuO^-s (NCCO) 
, has also shown a second peak transition; here, however, Bq decreases with temperature near the superconducting 
^0 ' transition at T c w 23K with no apparent multicritical poin& The second peak phenomena is also pronounced in 
other layered systems such as NbSe ^ 9 ^ 10 and in Pb/Ge multilayers^. Recent decoration datai 2 * on NbSe2 has shown 
that the topology of the vortex structure is weakly affected by crossing part of the second peak line. Hence the nature 
of the phase at B > Bo is not well established. 

The Josephson plasma resonance is a probe of the Josephson couplingi^*i4 and can be used to probe the various 
phase transitions. Recent data on BSCCO has indeed shown a significant reduction in the resonance frequency at the 
-j->; . second peak transitioni£*i&. 

In a remarkable experiment Fuchs et al^£ have shown that the phase diagram of BSSCO is much more elaborate. 
They show that the spatial distribution of an external current exhibits a transition from bulk pinning to surface 

1 ■ pinning of vortices with most of the current flowing at the sample edges. This depinning line crosses the multicritical 

point and its temperature is almost B independent at B < Bq. The depinning transition correlates with anomalies 
in vibrating reed experiments^ and in magnetization^. Thus there are four transition lines which emanate from the 
multicritical point at Bq, Tq: The first order line, the second peak line and depinning lines for both B < Bq and 
B > Bq. The common intersection of these four transition lines was also seen in data of the c axis Josephson critical 
current. This critical current decreases significantly above the second peak line (in contrast with the critical current 
^> , parallel to the layers) and also decreases in the depinned regimes. 

The notion of vortex matter in the presence of disorder has emerged as a fundamental problem of elastic manifolds 
in a random mediaSi. This has motivated an extensive theoretical effort towards understanding the field-temperature 
(B — T) phase diagram in presence of disorder. Impurity disorder does not allow long range translational order of the 
flux lattice and finite domains are expected 2 * 2 ^ At low temperatures and fields the system is a Bragg glas o 23 i 24 , i.e. 
the lattice is dislocation free, at long scales the displacement correlations decay as a power law and Bragg peaks are 
expected. The impurity induced domains are essential for the description of both equilibrium, e.g. thermodynamic 
phase transitions and non-equilibrium, e.g. critical current phenomena. Melting, e.g., is expected to occur by thermal 
or disorder induced dislocations, as indeed demonstrated for fields parallel to the layerS" 2 ^* 2 ^. Numerical simulations 
on related XY models have also shown disorder induced melting2LS2i2S. 

The flux lattice can undergo a transition which is unique to layered superconductors, i.e. a decoupling transition^S^i. 
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In this transition the Josephson coupling between layers vanishes while the lattice is maintained by the electro-magnetic 
coupling between layers. A disorder induced decoupling was also proposed as a crossover phenomena^. Decoupling 
in presence of columnar defects was also studied^ 3 -, showing enhancement of the coupled phase. 

It has been shown that decoupling coalesces with a defect unbinding transition- 4 *^ which has analogs in isotropic 
systems—. The resulting vacancies and interstitials lead to a reduction in the elastic tilt modulus^, consistent with 
the decoupling scenario as described below. It is possible then that a decoupling-defect transition accounts for the 
peak phenomena in all type II superconductors. The analysis below is, however, presented for layered anisotropic 
systems where quantitative predictions can be made. Vacancies and interstitials are neglected; their role is dicussed 
in the concluding section of the preceding companion article 38 . 

In the present work we expand our previous work^iMS and study effects of disorder at temperatures below the 
melting temperature T m by employing replica symmetry breaking (RSB) methods. The RSB methods are accurate 
when couplings of the nonlinear terms are weak. E.g. in the pure case they reproduce the RG result at weak 
Josephson coupling—; in the related problem of vacancies and interstitials it was shown that RSB accurately locates a 
disorder induced transition 3 -^. In the present problem weak coupling corresponds to weak Josephson coupling and weak 
disorder. Weak disorder can be stated as a condition on the size of domains Rbg being larger than the renormalized 
penetration length in the c direction. This condition is examined in section IV and the RSB actually detects this 
by producing a stronger singularity (appendices A-C). Furthermore, RSB as a variational method can identify order 
parameters and determine the form of the phase diagram. The critical behavior near the transition, however, is not 
expected to be accurate. 

The most interesting finding in this work is that of a glass order parameter which we term as Josephson glass (JG), 
as it is due to disorder induced on the Josephson coupling. The JG order is expected to lead to stronger pinning, 
hence the line where JG vanishes is associated with a depinning line. We find that the JG and decoupling lines 
cross and lead to four distinct phases which meet at one point in the B — T phase diagram, remarkably close to the 
experimental phase diagram 17,2 ^. This paper follows a companion oneS where the decoupling transition is studied in 
the pure system by second order renormalization group (RG). 

The full problem addressed here involves the following set of nonlinearities: (i) Josephson coupling which involves 
both pancake displacements and a nonsingular phase, (ii) A disordered Josephson coupling which leads to the JG 
order, (iii) A nonlinear coupling of disorder to the displacement pattern, leading to the well studied Bragg glass 
(BG) 23,24 . After presenting the model in section II, we study in section III a simplified version of the full problem in 
which the nonsingular phase is neglected and also the disorder coupling is linearized, corresponding to scales within 
finite domains. These approximations lead to an unphysical divergence of an integral I(z) where z is the renormalized 
Josephson coupling, i.e. z — > at decoupling. In section III we assume that I(z) is convergent and behaves as 
~ lnz, an assumption that is justified in appendices A, B and C. In appendix A we extend section III to solve the 
combined BG/JG system, though the nonsingular phase is neglected. In appendix B the BG system including the 
non-singular phase is solved, but JG is neglected, as relevant to thermal decoupling. In both appendices A and B 
we find an additional In 2 z term which signals a divergence of disorder effects in a regime close to decoupling. In 
appendix C we study JG with the nonsingular phase, but disorder is linearized. It is shown that I(z) converges even 
in this situation, while an additional \/\fz term is generated. In section IV we present a dimensional derivation of 
domain sizes which correctly reproduces the pinning and BG lengths. Near decoupling there is a regime of nonlinear 
elasticity with an apparent jump of the tilt modulus C44 and the critical current. This anharmonic regime coincides 
with the onset of the In 2 z term in appendices A, B. In section V the Josephson plasma frequency is studied, being 
an efficient probe for identifying the various phases. In section VI we discuss available data on the second peak and 
depinning transitions. We propose that decoupling accounts for the main features of the second peak transition while 
the depinning transitions correspond to the onset of JG order. 



Consider a flux lattice with an equilibrium position of the Z-th flux line at vectors R; of a regular two-dimensional 
lattice. The flux line is composed of a sequence of singular points, or "pancake" vortices, whose positions at the 
71-th layer can fluctuate to R; + u™ . Of particular interest is the transverse part of u" with the Fourier transform 
ut(<\, fc), where q, k are wavevectors parallel and perpendicular to the layers, respectively. The elastic energy due to 
the electromagnetic coupling has the form 



II. THE MODEL 




(1) 
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where the flux line density is 1/a 2 , d is the spacing between layers, q is within the Brillouin zone [of area (27r/a) 2 ], 
|fc| < n/d and k z — (2/d) sin(/cd/2). The shear and tilt moduli are given (for a> d) by 41 i 42 i 43 

c° 6 = r/(16da 2 ) 

c%(k) = [T/{Ma 2 \l b kl)]\n{l + a 2 kl/Av) (2) 

where r = (j^d / '(Att 2 A 2 b ) sets the energy scale and X a b is the magnetic penetration length parallel to the layers; 
t w 10 3 — 10 4 K for YBCO or BSCCO parameters 1 -. Note the strong dispersion of c1 4 (k) so that c® 4 (k) decreases by 
the large factor id/a) 2 when k varies from k < 1/a to 1/a < k < ixjd. 

The Josephson phase between the layers n and n + 1 at position r in both layers involves contributions from a 
nonsingular component and from singular vortex terms. The singular phase around a pancake vortex at position 
R; + u" is a(r — R; — u") where a(r) = arctan(y/x) with r = (x, y). We assume that all vortices belong to the flux 
lines, i.e. there are no free vacancies or intcrstitials. 

The Josephson phase involves the interlayer phase difference from the pancake singularities a(r — R/ — u") — a(r — 
R - uj l+1 ), which after expansion in u™ becomes (Eq. 19 of the companion article 3 ^) 



bn(r) = 2md f pM e _ n . r _ iknd(eikd _ 1)UT (q,fc) 
Jbz (271") q 



(3) 



We consider first a simplified model which neglects the nonsingular part of the Josephson phase. The nonsingular 
phase is essential for evaluating displacement fluctuations (section IV), however for the purpose of the phase transitions 
under study it can be neglected (justified by Appendices A,B). We have then an effective Hamiltonian for wavevectors 
M < Qo> (Eq. (23) of the companion article 3 ^) 

^i e /T=I^ C (< Z ,fc)g 2 |6(q,fc)| 2 ~|^^ / d 2 rcos6„(r) (4) 

where Ej is the interlayer Josephson coupling energy per unit area and 

^)-J^i^) + y^c(k) + 4 z (5) 
The last equality defines c(k) and c', i.e. 



c(k) 



to 2 ln(l + a 2 k 2 /4ir) 



32ir 2 Td 3 \ a 
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b 



Mn 2 Td 3 



(6) 



Since Va ~ 1/r decays slowly, even if u™ are small the contribution of many vortices which move in phase (q — > 0) 
leads to a divergent response of b n (r), i.e the 1/q factor in Eq. J2J. This leads to a decoupling transition^i* 3 ^* 3 ^, which 
at weak Ej is (Eqs. (27, 40) of the companion article^) 

_ 4a 4 f dk x _ to 2 log(a/d) 

d d 2 { J cUkY ~ 4^A 2 b ■ [ '> 

We note that melting and related dislocations have been neglected. An estimate of T m by the Lindemann criterion 
yields 2 ^ T m w r, hence our description near is limited to to a < A a (,. Melting in the absence of Josephson 
coupling was in fact studied 4 ^, showing that T m is between r/8 and the two-dimensional melting temperature of 
« 0.004r, approaching the latter at high fields a -C A a &. At intermediate fields the present description is then valid 
at a < OAXat- However, for disorder induced melting we estimate (see the discussion section VI) that for BSCCO 
parameters the decoupling field is below the melting field if a > 0.14A a b, consistent with the low temperature second 
peak field value. 

We proceed now to study the disorder term. A second assumption of the simplified version is that of linearized 
disorder, i.e. small fluctuations |up| <C a. Consider a short range pinning potential Up\ n (r) with the coupling 

H pm = fd 2 rJ2 U;\ n (r)p(r - R u?) (8) 
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where p(r) is a shape function for a vortex of size £o and the disorder has short range correlation 

(U; m (v)U;l n (v')) = \U5 n>n ,8{v r') (9) 

Expanding Eq. JBJ to first order in up and averaging £/™ in (r) by the replica metho d 2 ?' 4 ? leads to a disorder term in 
the free energy 

^2/T = |§££ur.ur' /3 do) 

n.l a. (3 

where a, f3 are replica indices. The average involves 

d lP {r)d lP (r)d 2 r = pS l3 (11) 
with p of order 1 . 

The replicated Hamiltonian of the simplified version, keeping only transverse displacements, is therefore 



n(1) / T = \ E [c(q,k)q 2 S a , -s Q yb a ( q ,k)b^(q,k) 

q,fc;a,/3 

"^E / d 2 r cosher) J d 2 r cos[6«(r) - &£(r)] (12) 



where 



2T 2 (27rd 2 ) 2 

It is found useful below to define a dimensionless disorder parameter s, 

r 2 o? \n 2 (a/d) 



UP a 2 d , . 

s a = ■ (13) 



(14) 



The inter-replica Josephson coupling, i.e. the £J„ term in Eq. I|12|) . is generated from the i?j term in second order 
renormalization group (RG). It is essential to keep the E v term from the start since it couples different replica indices 
and leads to distinct physics by RSB, as shown below. Physically, the E v term originates from random displacements 
of pancake vortices due to intralayer impurities. The pancake vortices are then not one on top of the other, resulting 
in random segments of Josephson vortices, i.e. vortices parallel to the layers. The latter represents random Josephson 
phases, whose replica average leads to the E v term in Eq. (|12f> . 

We proceed now to present the full model, which extends Eq. (|12l) to include the nonsingular phase as well as 
nonlinear disorder. The Josephson phase involves a nonsingular phase n (r) in addition to the pancake fluctuations 
via b n (r). The Hamiltonian of the pure system is then (Eq. 21 of the companion article^) 

H pure {b, 6}/T = i G7 X (q, k)\0(q, k)\ 2 + ± ]T c(q, fc)g 2 |6(q, fc)| 2 

q,k q,/c 

"fr£ fd 2 rcos[9 n (r)+b n (r)} (15) 

where 

Gf(q,k) = ^^(X- b 2 + kl). (16) 

Consider next the general form of the disorder coupling^!. Using the relation ^ 5 2 {p — R;) = J2i e l Qi 'P where Q ; 
are reciprocal lattice vectors, the disorder coupling JSJ) becomes 



it'pi'n 



J d 2 r E^inW / ^VV P u»(p)] E elQ ' P (1 ? ) 
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For |Q ; | < l/£ we can replace p(r) by £o<5 2 (r) so that 

^ n I 

The coupling to long wavelength modes via Vu"(r) is irrelevant in 3D so that the replica average of Tt P i n becomes 

n dis /T = ^ J d2r C0S tQ • ( U " Q ( r ) - ( r ))l ( 19 ) 

Q,a,0,n 

with go = US^o /T 2 a 2 . To relate this form to the linearized one (|10f) we expand in u n,a (r) and use 5Zq Q 2 ~ 
f- J 1 ^ Q 2 dQ 2 w a 2 /87r£g so that (|10|l is obtained if p w 1/27T. The coupling g can then be written as 

■ 9 ° " TO S 4T 2 A ab 4 ■ (20) 

We are interested here in BG effects on the q — ► singularity associated with the decoupling transition, i.e. the long 
range properties of the BG. The BG domain size is defined by the scale R where the displacement correlation starts 
to diverge as lnr. It is reasonable to expect that this scale is determined by the shortest Q, as indeed shown for a 
system with regular elasticity 2 ^, i.e. far from decoupling. We consider then the disorder term with just the shortest 
reciprocal wavevectors |Q| « 2ir/a (e.g. six wavevectors in the hexagonal lattice). The full Hamiltonian is then 

n / T = J2^ n p^{b a ,0 a }-^ Yl /^cosK(r)-^(r) + C(r)-^(r)] 

~ /rf 2 ^cos[Q- (u™' Q (r)-u™^(r))]. (21) 

^ -La J 

We note finally that a similar two-dimensional (2D) model has been studied by RSB and RG methods 45 ' 46 . As 
shown in the next section, finite values of k dominate the phase transitions, so that a certain k averages of the 
coefficients in Eq. (I12|) lead to a 2D problem with the same q singularities as in (|12|l . Indeed the RSB solution below 
has the same structure as the 2D case^ 5 - with a temperature parameter t — T/T® and a disorder parameter s (Eq. 
I14J) . In view of this similarity, it is useful to quote the RG equations of the 2D model 4 ^ in terms of u — £ 2 Ej/T and 

v = esjT, 

du = [2u{\ - t - s) - 2~f'uvt\d ln£ 

dv = [2v(l - 2t) + ± 7 'su 2 - 2itv 2 ]d ln£ 

dt = -2 7 " 2 (t + s)t 2 u 2 d\n£ 

d(s/t 2 ) = 167"Wdln£ (22) 

where the initial value of the scale £ is a and 7 ', 7" are numbers of order 1. We quote these results so that the necessity 
of the E v term is shown more concretely. Indeed E v is generated by sE 2 while at t < ^ it is relevant on its own. 
Furthermore, the RG results will be used to qualitatively support and supplement the phase diagram, as derived by 
RSB in the next section. 



III. PHASE DIAGRAM 



In this section we consider the simplified version, Eq. (|12|) . This assumes that displacements are within finite 
domains and Bragg glass effects are neglected; also the nonsingular phase is neglected here. Appendices A and B 
show that these assumptions are justified for the purpose of our phase diagram. The nonsingular phase is essential 
for evaluating displacement fluctuations, as studied in Appendix C. 

We proceed by using the RSB method^!. The RSB method proceeds by employing a variational free energy 
J~ var = J~o+ < H — Ho > with To the free energy corresponding to 



«o=i^ G" 1 /3 (q,A ; )6 Q (q,fc)6' 3 *(q,fc) 

q,fc;a,/3 



(23) 
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and G Qi( g(q, k) is determined by an extremum condition on F var . We define the following averages (...)o with respect 
to H , 



(cos6«(r)) = e-2 A " 

A a = G a:a (q, k) 



(cosK(r)-&f(r))o _^ 

E a ,f3 = 2^[G Q , Q (q,fc)~G Q , (q,fc)] 



(24) 



(25) 



q , k 



so that 



= hJ2 ^[lnGfo *>) + (G^fa, k) ~ c{q, k)q 2 I - S(] ^L)G{q, k)] 



q,fc 
E.j \ - 



Ey V ^ 



(26) 



where I a> p = 5 Q and L a p = 1. 

The variational equation SF var /SG a .p — yields 



G a]^ k ) = [c{q,k)q 2 + z\5 aJi - s a {q 2 /k 2 z ) - a a ^ 



a a,f3 



Ej I 4 



(27) 
(28) 
(29) 



where z is a renormalized Josephson coupling. In the replica limit with the number of replicas n — ■> the RSE 47 
method represents each matrix as a hierarchy of matrices, e.g. <y a ,f3 is represented by <r(u), with < u < 1 and a 
diagonal component a. We parameterize therefore G~^ by a and a(it), where 



a = c(g, fc)g 2 — s — + z — a 



a(u) 



1 i \ 
-so-^-a(u). 



The amount by which the replica symmetry is broken is measured by a glass order parameter 

r 

A[u) — ucr(u) — / cr(v) dv . 
Jo 

The inverse matrix G a ^ is represented by 6 and b(u), where^i (see also Appendix B of Ref. |4 



(30) 



(31) 



b = 



1 



-o(0) 



1 



dv A(v) 



b - b{u) 



v 2 a - (a) + A(v) 
1 dv 1 



u[a - (a) + A(u)] J u v 2 d-(a)+A(v) 



(32) 



and 



(a) = J a{v)dv = s — - (a) 

a-(a)+A(u) = c(q,k)q 2 + z + A(u) . 
B(u) can be written, using i|25|l and the inversion formula (|32|) . as 



(33) 



(34) 
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where 

c(k)q 2 + z + A(w) J 27r 47rc(fc) z + A(u) 



q,fc 
and 

„ f dk , c(k) 

c ^j^W)^MdY (36) 

Here c(<7, fc) of Eq. (|SJ) is replaced by c(fc) as denned in (J21 while the q 2 term in Eq. (J5J amounts to redefining the 
upper cutoff into q 2 = A\n{a/ d) / \ 2 ab , (considering k ~ 7r/e? as the dominant range of the following k integration) and 
A c = c(ir/d)q 2 . In the following a variable t is temperature in units of the pure system (Eq. UJ, i.e. 

_ T _ f dk 1 



Eq. (|2S1) is then 



g ( U ) = 2fln c . + Ci . (38) 
z + A(u) 



To find A(u) we note that Eq. I|29|) is equivalent to a(u) — (E v /Td) exp[— B(u)/2]. Differentiating this equation 
and using A'(u) = ua'{u) we obtain 



which by using (|38[) can be written as 



u du g (u) 



,1 LdA 

(--^)-r = 0. (40) 
u it du 



The solution of this equation is a one step function, i.e. A(u) jumps at u = 2t from zero to a constant value Ao at 
2t < u < 1. The solution is therefore nontrivial if t < 1/2. 

To complete the solution, the function B{u) from P4f is needed 

\B(u) = Ci + In — + (2t - 1) In Ac A u < 2t 
z z + A 

A r 

= C*i+2tln — 2t<u<l (41) 



which yields for cr(u) 



z + A 



z 

a(u) = (To = 1— o\ u <2t 

Z + An 



TcT A 

Finally, from Eq. 1)31(1 we have z + Aq = 2toi, hence, 



ai = |L(£ + ^0)» e -Cx 2t<«<l. (42) 



A c dT^' A 



b^"*)^ (43) 



A consistent weak coupling solution is indeed possible only at t < \. 

To find a second equation for z from l|28l) we need the first inversion formula in IMl'l) 



W ' ; c(% 2 + z (c(k)q 2 + z) 2 c{k)q 2 + z + A V ' 
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and after the q summation 



'■ = £ *<«■*> = / £1555 ^ + jay - ^ + ^5] + ft + £[/(•> + (45) 



q,fc 



where 



/" dq 2 dfc 1 

(Z) = 7 fc 2 c(fc) c(% 2 + z { ) 

and /'(z) = dl(z)/dz. For Ao ^ we have from Eq. 142|) (To = z/2t while for Ao = (possible at t > | as found 
below) we have Co ~ z2t "C 2, hence, with s defined in {3}, 

A a = ln(2et£; i ,/zd) + (so/87r 2 )[/(z) + z/ / (z)] A 7^ (47) 

= Ci + 2(t + s) In — + 2s A = . (48) 

Formally I(z) diverges at k = 0; this divergence can be traced back to our assumption that the cos[Q • (u"' Q — u™'' 3 )] 
term is expanded into the sq term in Eq. (|12|l . Retaining this cosine leads to domains of correlated u™. In Appendix 
A the joint BG-JG solution is found and is shown to remove the k — > divergence. A combined BG with non-singular 
phase solution is also shown in Appendix B to remove this divergence. The presence of BG, however, produces a 
term ~ In 2 z in a regime near decoupling (Eq. IA23I in the regime of Eq. IA27(I . This "anharmonic" regime is studied 
further in section IV. Excluding this anharmonic regime, the dominant part of I(z) is 

r /d dk f d <l 2 ird , A c 

The I(z) term in Eq. (|45|l can then be written as 



Therefore, the renormalized Josephson coupling of Eq. (j28(l is for Ao 7^ 0, using Eq. 1)47(1 . 

~A~ c ~ e ^TtdEvAj 1 ' 23 ' ^ 

Note that E v is generated from Ej by RQi^, i.e. E v ~ Ej initially; however, E v is RG relevant at t < \ even in 
1st order RG (Eq. I22|) . hence we consider E v and Ej as comparable so that E 2 j / (2TtdE v A c ) <C 1. Hence a consistent 
weak coupling z/A c <C 1 solution is possible only for s < h. Thus s — | marks a disorder induced decoupling with 
z = at s > |. 

Comparing Eqs. 1)4315 1|) shows that Ao vanishes at s = t (up to 0[ln(E j/E v )/ In E v ] term, small for Ej w E v 
dA c Td). Formally there is a solution with Ao < when s < t. However, the average distribution 47 of |6(q, fc)| 2 , which 
is ~ cxp[— |6(q, fc)| 2 /Go, ia (q, k)\, is acceptable only if G QjCe (q, k) > 0; this is therefore a thermodynamic stability 
criterion. Note in particular from Eq. (|44|l 

12 1 9f 

G aa (q = o,k) = -( —r-). (52) 

2t z z + Ao 

When s < t and Ao < the power dependence in Eq. I|43|) implies that z + Aq <C z (unless too close to t = s, i.e. 
s — t ~ 1/1 \nE v \) and therefore G aa (q — 0,k) < 0. This shows that only A > is acceptable. 

Thus the regime where both z, Aq are finite is limited to s < i, t < s; we term this regime the coupled Josephson 
Glass (JG) phase. The "coupled" notation means that the renormalized Josephson coupling is finite, i.e. z ^ 0. The 
glass parameter vanishes (continuously) at t — s while the Josephson coupling vanishes (with an apparent discontinuity 
- see section IV ) at s = h (see Fig. 1). For s > i and t < \ the solution is z = while A 7^ satisfies Eq. (JI3J, i.e. 
it is a decoupled JG phase. Recall that the JG order parameter Ao is due to E v which is initially generated by Ej. 
In fact, the RG of Eq. (1221) shows (see a similar effect in Fig. 3 of the companion article^ for the pure system) that 
Ej first increases (scaling from £ to l/q u ), generating the E v term, and only at scales beyond l/q u Ej decreases to 
zero. It is remarkable then that Ej is renormalized to zero while the JG order survives, much like the smile of the 
Cheshire cat. 
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Finally, for Ao = a replica symmetric solution is valid at s < t, which upon using Eqs. (|28l I48fl becomes 

Thus s + i = lfors<i defines a " thermal" decoupling transition. 

The interpretation of the phase diagram needs to be supplemented by a few observations from an RG analysis. The 
RSB results above coincide with those of a 2D model where the parameters t, s of the 3d system, as suitable k averages 
(Eqs. I37I5U|) . correspond to Hamiltonian parameters of the 2D system 4 ^. With this correspondence in mind, we infer 
next some qualitative modifications by using the 2D RG equation o 45 ! 46 , Eq. (|22|l . Note first that in a coupled phase 
z is RG relevant and therefore E v , which is generated to order z 2 , is finite too, hence a weak glass phase is expected 
also in the regime s < t < 1 — s; this weak glass order is not captured by the RSB solution. The line t — s for s < i 
can therefore be either a 1st order transition or a crossover line. RG suggests (Eq. I22|) this crossover line at t = \: at 
t < \ RG yields E v which is largely independent of z, hence a strong JG order, while at t > h RG generates E v ~ z 2 
with a weak JG order. The stability of the RSB solution shows that in fact this line, which is either 1st order or a 
crossover, is at t — s. 

TheRG, shows also a disorder induced decoupling, since Eq. \\2.'2\\ has a fixed point with u* = and v* = (l — 2t)/j't, 
stable at t < 5 and strong disorder. Note that for this solution s ~ ln£ increases with scale £, hence the correlator 
T(r) = (cos 6^(r) cos 6"(0)) which by RSB decays as r~ 2 ~ As is actually decaying faster as lnTfr) ~ — In 2 r. Explicit 
solution of the 2D RG equations^ found indeed a phase diagram very similar to that in Ref. |45| or in Fig. 1. 

The phase diagram, shown in Fig. 1, has three phase transition lines and a line which is either a 1st order or 
a crossover line. All these lines meet at a multicritical point s = t = 1/2. We interpret the transition where Ao 
vanishes as a depinning transition, i.e. the JG order parameter provides an additional pinning to that from the Bragg 
glass. The phase diagram has then a decoupling line which crosses a depinning line at the multicritical point. The 
decoupling line has a disorder driven section, s = \ aXt< i. 
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FIG. 1: Phase diagram. Full lines are decoupling lines where the Josephson coupling vanishes. The upper dashed line is a 
depinning transition where the Josephson glass parameter vanishes; the lower dashed line is either a 1st order line or a crossover 
into a weaker JG phase, i.e. weaker pinning. 
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FIG. 2: Phase diagram in terms of field and temperature. Full lines are decoupling lines [B = Bo and B = 2BoTo/(T + To)] 
where the Josephson coupling vanishes. The upper dashed line is a depinning transition (T = BqTq/B) where the Josephson 
glass parameter vanishes; the lower dashed line (T = To) is either a 1st order transition or a crossover into weaker pinning. 



The phase diagram in terms of field and temperature is derived by defining Bq,Tq as the field and temperature 
value of the multicritical point and is shown in Fig. 2. Bo is determined by the disorder strength via s = o while 
T = ±T d °(a = y/B /<j> ) (Eq. 0). Hence s = B/2B Q and t = TB/T B Q , up to InB terms. Since s increases with B 
the s = i line defines a decoupling transition from a coupled JG phase at low B to a decoupled JG phase at high 
fields. 

The coupled JG phase at B < B goes through either a 1st order or a crossover line at t = s, i.e. at T — T (up 
to InB factors). Therefore at T > Tq the glass parameter Ao is significantly reduced implying depinning, a change 
from strong to weak pinning. The decoupled JG phase undergoes a depinning transition into a decoupled phase at 
T = BqTq/B. Note that all phases, even the high T decoupled one, are Bragg glass phases of the flux lattice; in the 
decoupled phase the lattice is maintained by the interlayer electromagnetic coupling. 

The JG coupled phase at T > To undergoes a decoupling transition at t = 1 — s, i.e. B = 2i?o?o/(T + T ). This 
transition is continuous; the variational method of the pure system has been formally extended to higher J/T and 
found to be of first order—. As shown in the companion article^, the transition remains 2nd order when proper 2nd 
order RG is employed. Disorder, however, leads to an apparent discontinuity near decoupling, as discussed in the 
next section. 



IV. DOMAIN SIZES 



In this section we estimate various domain sizes and evaluate displacement fluctuations which identify these sizes. 
Remarkably, the expressions for the domain sizes are confirmed (up to numerical prefactors) by BG solutions (Ap- 
pendices A-C). The nonsingular phase, which was irrelevant for the purpose of the phase diagram in section III, is 
essential now. 

To appreciate the effect of the nonsingular phase 9, we briefly review the derivation of the transverse tilt modulus 
C44 of a pure flux lattice*. The Josephson phase involves the contribution of pancake fluctuations via b n (r) as well as 
a nonsingular phase, with the Hamiltonian Eq. i|15|l . To identify C44 we expand the Josephson coupling to 2nd order 
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in b(q,k) = 6(q, k) + 0(q, k), 



H pure /T= i^{[G7 1 (< Z ,fc) + ^] 

q,k 



6(q, k) - 6(q, k) - ; E 



Gf(q,k) + 



Td 



G f (q fc)% ^ fc)|2 + fc)|&( ^ fc)|2 + Q{¥ ^ k))} (M) 

The first term decouples from 6(q, k) and with |6(q, k)\ 2 = (2ird 2 ) 2 k 2 z \u tr (q, k)\ 2 /q 2 (Eq. 0) we identify 41 i 42 i 43 

cuiq , k) = cUk) + g 1 + Agg2 ' +A ^ + 1*V4*8> (55) 

where X 2 — A 2 b r/ (4ird 2 Ej); the last term is from reducing high momenta of the 2nd term of l|54|) into the 1st Brilluin 
zone. 

The second term of C44(g, k) is peculiar: at q ^ it vanishes when Ej vanishes and A c — > oo, as it should. However, 
at q — this term seems to survive even if A c — > oo. The origin of this peculiarity is that the harmonic expansion 
of the Josephson cosine term which identifies C44 fails 4 ^ when both q, 1/A C — > 0. The shift in the 1st term of Eq. 
(|54|l identifies an expansion parameter— with terms ~ q 2 k 2 \uT(q_, k)\ 2 /[q 2 + A,7 2 (l + X 2 b k 2 )} 2 , which diverge when 
both q, 1/A C — > and the expansion becomes invalid. In fact, the nonlinear cosine term replaces Ej /Td by z or A c is 
replaced by a renormalized 

Af = ^\ 2 ab r/{^Td3z) (56) 

which diverges at decoupling. Hence usual elasticity at q, 1/A,? — > near decoupling is ill defined. 

The Bragg glass domain size Rbg (parallel to the layers) sets a scale for the relevant q values. When Rbg > A,? the 
tilt modulus is large, containing the B 2 /4ir term of Eq. 1)55(1 . However, as decoupling at the field Bq is approached A^ 
diverges so that when Rbg < Af Eq. 155fl fails to describe C44 on the scale of q « 1/ Rbg- This defines an anharmonic 
crossover regime where usual elasticity cannot be used to derive Bragg glass properties. Finally, at B > Bg elasticity 
is restored and C44 is reduced to the first term in Eq. Q55JI. The main interest is in the regime of strong fields, i.e. 
a < 2A a f, where T <^ r is below melting. Thus at B < B and for sufficiently large domains the second term in Eq. 
(I55|l dominates and C44 = B 2 /An while at B > Bq only the magnetic coupling survives C44 = c^fc) which at ka 1 
becomes t / (Z2iT\ 2 ab d). Hence there is an apparent discontinuity, 

c 44 = 7rA 2 fc r/da 4 Af < R BG (57) 

= T/(32Tr\ 2 ab d) Af = 00 (58) 

Thus C44 is reduced within the anharmonic regime by the small factor 

e = a 4 /(32^ 2 A 4 h ) . (59) 

The apparent discontinuity in C44 affects also the domain sizes which can be estimated by a dimensional 
arguments 2 * 2 ^. Consider the tilt C44 and shear cqq terms of the elasticity Hamiltonian for the displacement u(r) 
and its transverse component uy(r). Rescaling parallel and perpendicular lengths yields an isotropic formSiiS4, which 
together with the pinning energy 1)18(1 yield (ignoring elasticity of longitudinal displacements) 

U= /"dM|c4i 3 c^[V U T(r)] 2 -(eoV«^)^n(r)^cosQ-[r-u(r)]} (60) 

J Q 

where the disorder coupling to Vur(r) is neglected. To estimate the energy gain from disorder we consider the 
overlap of the disorder energy between two configurations u(r) and u'(r) which are solutions for two realizations of 
the random potential^ 1 .; this overlap is a measure of the energy variance in configuration space. The r integration 
leads to a single Q sum so that the variance is ~ ^2q cos Q ' [ u ( r ) — u '( r )]- Each of u(r) and u'(r) has fluctuations 
(u 2 ) ps {uj) in a domain of size R' so that the Q sum is cutoff by Q < (w 2 ^) -1 / 2 . Below this cutoff the cosine can be 
expanded and summed so that averaging Eq. I(6U|) yields 

(H)/R« = \c^cli\u 2 T )R'- 2 - U^ 2 e /[aMu 2 T )R' 3 } 1/2 ■ (61) 
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Minimizing with respect to R' yields R 1 ~ (u^) 3 , i-e. the Flory exponent. This exponent is not exact; the more 
accurate statement, shown within the BG solution? 3 , is that the disorder averaged correlation B(R') = ([ttr(R') — 
mt(0)] 2 ) ~ i?' 1 / 3 is a quantitatively correct description in the range between the pinning length R p where B(R P ) = £ 2 
and Rbg where B(Rbg) = fl2 - The fluctuations (uj>) on scale i?' in the dimensional argument correspond then to 
B(R'). 




FIG. 3: Bragg glass domain size Rbg parallel to the layers and the renormalized London length perpendicular to the layers 
A,?; the latter diverges at the decoupling field Bo- Rbg can be found from elasticity for B < Bo only if Rbg > A,?; otherwise, 
as in the hatched region, the elastic tilt modulus is ill defined. 



The domain size parallel to the layers is, from minimizing Eq. (|61|l . (up to \n(a/d) and a numerical prefactor) 

R+ « (A Qb /a) 5 (4) 3 /( S £ 4 d) \?<R+ 

R- « (\ ab /a) 3 (u 2 T ) 3 /{4ir S $d) A* = 00. (62) 

The pinning length R = R p is given by Eq. (B2J| with {u 2 T ) w Co- The condition Af < i?+ is not valid for BSCCO 
parameters; to allow for large pinning domains one needs either a <C X a b or to allow for domains with a somewhat 
larger fluctuations in {uj); the latter increases R p very rapidly since it increases with the 6-th power of ut- The 
critical current can now be estimated2ii2£ by balancing the Lorenz force j c BR 3 /c with the pinning force (H)/£o 
(evaluated at the minimum of Eq. (|61[) h leading to j c ~ I/C44. Increasing the field within the anharmonic regime 
decreases C44 by the factor e so that j c is enhanced by a 1/e factor which is significant when a < X a b- 

A second length scale R = Rbg is identified by Eq. l|rj2^ with the fluctuations (itf,) w a 2 . The proper definition of 
Rbg is the scale for the onset of the lnr form for the displacement correlation function, as inferred in Eq. l|A26jl or 
It is remarkable that Eq. (|62|l gives the correct form for for Rbg, up to a numerical prefactor, i.e. Eqs. I|A26I 
Eq. (|62|l shows that Rbg is reduced by e 1 / 2 through the anharmonic regime. The onset of the anharmonic 
regime is at Rbg ~ i- e - 

X? „ 10- 3 ^ (63) 

with a numerical prefactor from the BG solution (Eqs. IA26I lB8|) . For BSCCO or YBCO parameters at s w i this 
reduces to A^/A a & ~ 10 5 , i.e. the initial anisotropy of X c /X a b = 10 — 100 has to increase to w 10 5 . Since z is 
exponentially renormalized (Eqs. 1511 153|) this anharmonic range may be observable. 

Fig. 3 illustrates the lengths Rbg and A,?, demonstrating the anharmonic regime within which Rbg has a significant 
drop and correspondingly j c has an apparent jump. Note that even in the decoupled phase (B > Bq) Rbg is large 
for typical type II superconductors, Rbg ~ ^ab - 3 1 '(^ 7rs ^o < ^) ^ a , consistent with a decoupling transition within the 
Bragg glass phase, i.e. below a melting transition. 
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The solution of section III can also be extended to include the nonsingular phase. Since disorder is linearized, the 
pinning length R p can be determined, though the BG length cannot. Appendix C develops this solution and shows 
that in the coupled phase R p coincides with R + Eq. (with (uj) ~ £o)j u P to a numerical prefactor. 

The main result is then that the fluctuations in Ur{r) behave with an effective C44 which is large when q < 1/A^ 
(Eq. 1)57^1. i.e. for domain sizes Rbg > A^, while for z = C44 is reduced (Eq. (J5SJ). While the condition A,? < R p in 
Eq. (|62|) is not valid for BSCCO (the pinning domains are likely to be two dimensional) our results for the anharmonic 
regime itself in terms of the much larger Rbg are valid. The existence of a narrow anharmonic regime leads to an 
apparent jump in C44 which possibly affects the critical current. 

In the anharmonic region below decoupling (see Fig. 3) where Rbg < A,? a more complete form [e.g. Eq. l|C8(l ] 
is required to interpolate between the limiting forms of c 44 . However, a method relying on an effective harmonic 
theory, such as RSB, is suspect within the anharmonic regime, since the system has no effective elastic constants. 
Furthermore, RSB signals this deficiency by producing a ln 2 z term, precisely in the the anharmonic regime found 
here, as shown in Appendices A and B. 



V. JOSEPHSON PLASMA RESONANCE 

Josephson plasma resonance provides extremely useful data for identifying phases of vortex matteniiiiiSiiS. In 
particular a jump in the resonance frequency u) p i has shown^i^ that the Josephson coupling is strongly modified 
at the second peak transition. In this section we derive u> p i in the ordered phases and also consider the fluctuation 
contribution in the disordered phase. The Josephson plasma frequency is given byi^ (see also the companion article^ 
section V) 

<4l= eQh2 (cosfe) (64) 

where eo is a dielectric constant. The task is then to evaluate the thermodynamic average (cos b). 

Consider first the ordered phases where at least one of z and A is finite. We start by evaluating T var of Eq. 
for a general one step RSB, recover the solution of section III, and then identify (cos b). This derivation is needed so 
that the free energy itself can be evaluated, and from the latter (cos b) is inferred. The self mass term a a ,p of Eq. 
H27J1 is written for a one step solution in the form 

a = <7 L + (01 - <tq)C - [aon + (ai — <J )m]I (65) 

where £ a ,/3 = 1 and C has 1 elements in blocks of size m x m sitting consecutively along the diagonal, and elements 
otherwise. For mOwc identify (o\ — uo)m = Aq so that 

k 2 m 



G{q, k) = [(c(q, k)q 2 + z + A )I + (-sofj - - ZT^ 1 =cJ + /3L + jC. (66) 



It is straightforward to identify the coefficients of the inverse matrix 

a(q,k) = — — , a = > a(q, k) = 2tln— — — + d 

c{q, k)q 2 + z + A e-J z + A 

2 

SOTT + do A 2tan 

(3(q,k) = I' /3 = ^/3g,fc)=2sln-^-l +— ~ 
[c(q, kjq^ + z\ z z z 

j{q,k) = - — [ a (q,k) - — 1 ] , 7 = V j(q, k) = — In Z + A ° (67) 
m c(q, k)q z + z m z 

where the form of Eq. (|50() is used for I(z) in the 2nd line. The definition of a identifies o\ = 2E v e~ a /d and 
(To = 2£ , 1 ,e _a_7 /<i. We follow a similar algebra in section IV of Ref. Ha to evaluate the free energy density per replica 
as 

f(m, z, A ) = /o + (1 - -tA + (1 + S -)tz ^(2t - m) e — 7 + ^(1 - m)e a EL e ~^^ (68) 
m t 1 a la la 

where /o is m, z and Ao independent. Minimizing /(to, z, Aq) yields m = 2t and Eqs. 142l43l51f) for uq, <j\, z + Ao 
and z. A replica symmetric solution is also possible with Ao = leading to Eq. I|53|l . The free energy at minimum is 

fmin = f Q + (t-l + — )(A + Z) + (S- \)z . (69) 



14 



The Hamiltonian Eq. l|12f) shows that (cos&) = —Td(df /dEj). As discussed below Eq. (|53|l E v is generated from 
Ej in 2nd order RG so that E v ~ E 2 initially, while E v is RG relevant at t < |, so that its value which is to be used 
by the variational scheme is more weakly Ej dependent. We assume then E v ~ Ej with < k < 2. Hence in the 
Ao 7^ phases 

<9(z + A ) k(z + A ) 



dEj Ej{\- 2t) 

dz (2 - k) 

~dE~j ~ (1 - 2s)Ej 

(cosb) = -«(1 - 2t) Z + A ° + (1 - \k)-?— (70) 

Zbare Zbare 



where Zb are = Ej/Td is the bare value of z. For the A = phase 

dz z 



dEj (1-t- s)E. 



j 



z 



(cos b) = (71) 

Zbare 

so that at T = s = the order is maximal, (cos b) = 1. 

These results show that the JG order produces a negative contribution to (cos b) so that when crossing a depinning 
line (cos6) is enhanced by the ~ n terms in Eq. (|70|) . Since Ao is continuous, the jump at depinning is fc(|— 2t)z/zbare- 
As discusses in section III, the depinning in the lower part of Fig. 1 is not a strict phase transition, but rather a 
crossover line, hence we expect a smeared jump of (cos b). An observation of a (cos b) enhancement when crossing the 
lower depinning line at T m Tq {B < Bo) would be a clear signature that depinning relates to JG order. The actual 
enhancement depends on k, for which we do not have a precise derivation. 

Near the decoupling transitions, the presence of anharmonic regimes, shown in section IV, lead to an apparent 
jump in (cos b). This jump relates to the z terms in (|TC)|) and also depends on the fluctuation contribution which is 
considered next. 

We proceed to evaluate fluctuation contribution when (cos b) is small. As shown by Koshele\ii^ the local (cosb n (r)) 
is finite even at high temperatures, e.g. above the decoupling transition. The high temperature expansion, while 
formally ill defined, does reproduce the RG results for (cos b), as shown in section III of the companion article 38 . The 
high temperature expansion yields 



(cos6„(r)) = (Ej/2T) J d 2 r exp[-A(r)] 



A(r) = £(l-cosq.r)(r(q,fc)| 2 >. (72) 

For r 2 > 1 /q 2 we can use the form (|45H with z replaced by a cutoff c(k)/r 2 while for r < l/q u we expand 1 — cos q • r — > 
\q 2 r 2 : hence 

A(r) = 4(s + t) Hq u r) r > l/q u 

= l(s + t) q y r<l/q u . (73) 

The two regimes in Eq. (|73|l give comparable results, though the r > l/q u is larger near the transition and 
reproduces the form of the RG result, as discussed in section III of the companion article 38 . The latter yields, in 
terms of the multicritical point coordinates (up to In B terms) , 

( cos b) pb ^ jA ° b ^ (74) 

[ 1 81n(o/d) T(BT + BT - B T ) ' 1 ' 

Well above decoupling at s + t h we obtain (cos6) ~ [BT(T + T )] _1 . A 1/BT dependence has been obtained 
by Koshele\iii with a weakly temperature dependent prefactor for an XY model, i.e. infinite X a b model. This result 
corresponds, in fact, to the melted, or liquid phase 38 . Data on BSCCO— has shown that (cos&„(r)) ~ £-o.8y-i m 
reasonable agreement with the 1/BT form. The present result shows that in the decoupled phase, below melting, 
u p i ~ [BT(T+To)] _1 , or the form l|74|) near decoupling. This distinct temperature dependence can be used to identify 
the decoupled phase. 

As decoupling is crossed, we expect a positive fluctuation term to compensate the negative contribution of the JG 
order. Thus the forms l|7UI71|l can be used for the jumps of (cos b) across depinning or decoupling, while (|72|l is valid 
in the high temperature or high field regime where (cos b) is small. 
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VI. DISCUSSION 

The present work exhibits the JG order parameter as well as the decoupling transition with disorder. We discuss 
now our proposal for each of the 4 transition lines emanating from the multi-critical point (Fig. 2) and compare with 
experimental data. 

Consider first the decoupling transition within the JG phase at B = Bo, T < Tq. We have shown that RSB methods 
are suspect within a narrow region near decoupling, where usual elasticity is ill defined (Fig. 3). RSB identifies this 
as a In 2 z divergence in A a which renormalizes z (Eq. I28J1 . This can be thought of as a disorder term s e f / ~ In z with 
a diverging s e //. The consequence is an apparent discontinuity, or even an intrinsic 1st order transition, driven by 
disorder. 

This decoupling transition is consistent with the main features of the second peak transition: (i) decoupling field 
being weakly T dependen t 1 ' 2 '?' 4 , (ii) decoupling field decreasing with impurity concentration^, (ii) an apparent jump 
in the critical currenlsiiSi&iiSaLifliii (iv) decrease in the c axis critical current and (v) a jump in the Josephson 
plasma resonanc o 15 ' 16 . The anharmonic region near decoupling leads to an apparent reduction in C44. The reduction 
in C44 and the resulting reduction in domain sizes account qualitatively for the enhanced j c . We do not attempt 
a quantitative fit; in fact, the measured magnetization changes (and inferred j c ) at the second peak decrease with 
temperature due to the strongly temperature dependent relaxation rates^, approaching the much smaller equilibrium 
magnetizations . 

The nature of the phase at fields above the second peak line has not been conclusively settled. This work proposes 
that it is a BG phase where the domain sizes have been reduced by ^fe — a 2 /(47r\/2A a 6 2 ). Experimentally, the smooth 
connection of the second peak with the 1st order line^ suggests that it is a single "order-disorder" line of common 
origin, e.g. a melting line. However, the presence of a depinning line that crosses the "order-disorder" line has been 
seen by numerous experiments 1 ii2*i2i2&. The crossing of this depinning line with the "order-disorder" line, separates 
the latter into a disorder driven second peak part within a pinned regime and into a thermally driven part in a 
depinncd, or more weakly pinned regime. This depinning line corresponds to the onset of a Josephson glass order, as 
suggested below. 

Consider next the decoupling line at T > Tq. This corresponds to the 1st order transition, which is considered as 
a melting linei^. However, neutron datai has shown a reentrant behavior in the 600 — 10 3 G range with positional 
correlations increasing with temperature. It is possible then that near the multi-critical point the 1st order line is a 
decoupling line. At higher temperatures decoupling then merges into a melting line. 

The 3rd transition line is a transition within the JG order at T — To, B < Bq into a weaker JG at T > Tq. A 
depinning line which is almost vertical at T w To was indeed observediii2ii2i2S. We note in particular the c axis 
critical current^ which shows a decrease on the high temperature side of the depinning line. The thermodynamic 
critical current is proportional to the renormalized z that changes from the weakly T dependent Eq. (51) at T < Tq 
to the strong exponential decrease with T in Eq. (53) at T > To , consistent with the data. At T > To we also expect 
a sharp enhancement of the Josephson plasma resonance, which is an additional tool for identifying the JG order 
parameter. 

The final 4th line is a depinning line at T = BqTq/B, B > Bq corresponding to a depinning line as observed in 
BSCCO by current distribution dataJ^, vibration reedi&, magnetization 1 ^ and c axis critical current data^S. This line 
is more difficult to detect by Josephson plasma resonance since its frequency varies continuously, with discontinuities 
in derivatives. In the decoupled phases (with or without JG order), where (cos&) is small, we expect the fluctuation 
form Eq. f71|). 

We have assumed throughout that our transition lines are well below melting. Thermal melting is discussed 
below Eq. Q while here we estimate the disorder induced melting field. We assume a Lindeman criterion such 
that the fluctuations in the decoupled phase on scale R~ = a are (u 2 ^) = cj^a 2 , with cl = 0.15 a conventional 
Lindeman number 21 . Using Eq. H62[) with a prefactor as identified by Eq. I|A26|) yields a melting field of B m w 
10~ 2 c 3 L y/2B B\\l b /d£,Q where B x = 0o/A 2 b . With BSCCO parameters the condition B rn > B is satisfied if B < 
50B\ , hence with the second peak field of B a w B x ~ 500G disorder induced melting is expected at a higher field. 

In conclusion, we have found a phase diagram which is remarkably close to the experimental on ei*2iiL2iiLi2i2S, having 
a multicritical point and providing a fundamental interpretation of both the second peak transition and the more 
recently observed depinning transitions. 

Note Added 

In a recent work [H. Beidenkopf, N. Avraham, Y. Myasoedov, H. Shtrikman, E. Zeldov, and T. Tamegai (un- 
published)] the depinning transitions were identified by relaxed magnetization data as equilibrium transitions. Both 
transitions at fields below and above the multicritical point were identified and suggested to be equilibrium glass 
transitions. 
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APPENDIX A: BRAGG AND JOSEPHSON GLASSES 



This section studies nonlinearities due to both disorder and Josephson coupling leading to two glass order parameters 
- the Josephson glass (JG) and the Bragg glass (BG); the non-singular phase is neglected. In particular an equivalent 
term to the integral I(z) (Eq. I46|) is identified and is shown to be convergent at k — ► 0. 

We consider the full Hamiltonian Eq. 121L which by neglecting the nonsingular phase becomes 



n/T=\ £( c '|! + c (%2)|&«( q ,A0| 2 -^£ |d 2 rcos6«(r) 

J2 J d 2 r coB[l£(r) - &£(r)] ~ % j ^cos[Q • (u"-(r) - u"^(r))] (Al) 



The average of the disorder term over the variational Hamiltonian (|23f) TCq yields 

(cos[Q • (u"<«(r) - u^(r))]) = exp{-^ ]T ^[G aa {q, k) - G a0 (q, k)]} . (A2) 

q,jfc 2 

We assume for simplicity a square lattice, Q = 2ir/a, otherwise a factor (aQ/27r) 2 is needed in the exponent; there 
are then 4 shortest Q terms in Eq. IjAlfl . The variational equation for G~^(g, fc), Eq. 1|27|) . has now an additional 

self energy term cr^g which allows for an additional RSB. Written as an equation for matrices in replica space, e.g. 

G, we have 

G-\q, k) = (c'fj + c(k)q 2 + z)i -a 2 - f^l • (A3) 



When (cti) q /3 = 1, i.e. no RSB, the previous form l|27|) is recovered. The variational T var <|26[) has now a term 
~ exp[— \B^1] (instead of the sq term) where 

B S = J E fi [ G «« (9. fc ) ~ G «P (iM- ( A4) 



The variation of this term identifies 



7 

while cr^g and B^g have the previous forms l|25l 129(1 . In the hierarchical scheme G -1 is represented by [5, a(u)] which 
are now given by 

A ? 

i y - /7 \ 2 r - - r 

a = c + c(k)q + z - cr 2 - (Ti 
q 2 

a(tt) = -o- 2 (u) - cti(u)— . (A6) 

The JG and BG order parameters which measure the degree of RSB are Ai(u), A2(u), respectively, where A, (it) = 
u<Ji{u) — Jq <Ji(v)dv, i — 1,2. Using the inversion l|32|l we can write 

i RM 9i(u) f 1 9i{v) -to i Av^ 

5-BiW = / — 5-«« i = l,2 (A7) 

u J u v 2 

where 

= ^E[ c V + c(fc)fc 2 + (z + A 2 ( M ))4 + A 1 ( M )]- 1 

S2(u) = Yyt> + c ( fc )« 2 + z + A ^ + M^Sr 1 ■ ( A8 ) 

q,fc z z 
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As in Eq. we find 

A>(u) d A>(u) 
u du g[{u) 



1,2. (A9) 



Consider first g2{u) which is dominated by k 3> q so that the c' term produces just the cutoff q u . The q integration 
then yields Eq. I|38l) with A c — > [c(k) + Ai(u) / k 2 ]q 2 in the logarithm. As above, we replace k by n/d in this logarithm 
since the k integral is dominated by k « ir/d due to the significant softening of c(k) near k = ir/d. Hence the form 
92 (u) ~ ln[z + A 2 (it)] is maintained and the solution, as in (|40|l is a one step function at u = 2t. 
To solve the equation for Ai (it) we simplify the form of c(k) as 

c(k) = c(0) = c_ k < - 

a 

'"OJil k> \ < A10 > 

This form captures the significant dispersion of c(k) with c(n/d) -C c(0) and allows analytic treatment of the potentially 
divergent k — > integrals. The /c > - integration range in gi(u) has an integrand 

[cV + A x (u) + c(vr/d)(4/d 2 ) + (z + A^u))*; 2 /,? 2 ]- 1 

so that c(ir / 'rf)(4/ 'd 2 ) 3> Ai (it) provides a cutoff on the g integration, i.e. <?i(it) acquires a term independent of Ai (it). 
The k < i integration has c(0)/a 2 3> Ai(u) so that after the k integration 

a 2 f q 2 do 

9i{u) = tt-^ / ; + const . (All) 

8rf J ^/( c > q 2 + A!(ii))(c_g 2 + z + A 2 (u)) 

Ai(u) varies between Ai(0) = and Ai(u c ) which depends on the disorder strength (see below); Ai(«) = Ai(u c ) 
is constant at u > u c , being a valid solution of (|A9|) . As the decoupling transition is approached and z — > the 
q integration in (|A11|) has distinct forms depending on the ratio of Ai (it) and c'z/c_. When Ai (11) < c'z/c- the 
dominant integration range is Ai(u)/c' < q 2 < z/c- and the result for the derivative is 

'S 1 (A 1 ) = ^lnA 1 a' = _^_. (A12) 



dA^ 1 " 1J y/i 1 87rdV 3 / 2 
Substituting in (|A9(I yields Ai (it) w Uy/z/a'ln 2 u and with l|A7jl we obtain [C2 = 51(11 = 0)] 

m it In it 



hBxiu) = C, + lnln?i + 0(-^-) (A13) 

In u 



so that CTi(it) ~ 1/ In u — > at 11 — > 0. When Ai(u) > c'z/c(0) the dominant integration range is z/c(0) < Ai(u)/d < 
q 2 so that 2 = can be taken and 

d , A n ct~ a 2 ? k -t a\ 

ffi(A0 = = = — . (A14) 



dA^ - VAT 8 ™'ci /2 

Substituting in (|A9|) yields Ai (it), so that in both regimes we have to leading order in u 

a / \ %-Jz . , . c' 

Ai(«) = ,/ 2 A 1 (tt)<z— 
a In u c- 

u 2 A , x d 



Aa 2 _ ly ' ' ~c_ 

Integrating (|A14|) and using for 2a_ A /c'z/c_ < it < it c , 

51 (it) = C 2 - u 



A 1 (u)>z—. (A15) 
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±£1(11) = Ca-ln u c (A16) 

llr 
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so that <Ji(u) ~ it in this range. We suspect that the solution at Ai(u) < zc'/c- is significantly modified by the 
non-singular phase (as indeed found in Appendix B). This is of no concern since anyway the effect of this range on 
the z equation vanishes (Eq. IA20l below). 

We finally consider the equation for z by using the inversion formula IMl'l) 



A a = / G(q, k) 



0- 2 (O) + <n(O)fr 



q,fc ' Jq,k c'|r + c(k)q 2 + z c'fy + c(k)q 2 + z 

l dv A 2 (v) + A 1 (v)£ 



o v 2 d£+c(k)q 2 +z + A 2 (v)+A 1 (v)£ 



(A17) 



Taking (72(0) ~ z from section III the (72(0) term yields a constant, independent of z. Note that without BG order, 
Ai(u) = 0, the one step solution for A 2 (u) reproduces the so terms in Eq. I|44fl. 

Consider first the range k < l/a which led to an apparent divergence in section III. For small v, where the v integral 
may diverge, we take A 2 (w) = so that 

Al = l ^iqlkKl/a ^cV + (c-q 2 + z)k 2 ~ d<f + (c_q 2 + z)k 2 + Al (u)ql ] (M8) 
Performing the k integral leads to a [c-q 2 + z]~ 3 ^ 2 factor, which amounts to a lower cutoff yJzJcZ, 

A ^= I t/2 I rf?[-V?+-VAiW+c'g 2 ] (A19) 
J 4ttcI / V J \fz~fcZ q 

For Ai(v) < zc' / C- one can expand in A\(v), which from l|A15|) yields a term 

/ ^T^^T^^O z^O. (A20) 
Jo v In v hi z 

For the v integration range where Ax(v) > zd/c_, which exists if Ai(u c ) > zd/c_, we have 

The second contribution to A a is from the range k > l/a where c{k)k 2 ~ constant provides a cutoff in the A a 
integrations, hence A\(v) can be neglected in the denominator, leading to 

A*= [ [ "I . .., / (A22) 



k>i/a Jq k 2 [c(k)q 2 + z} 2 

Identifying so = / $Ai(v) we obtain the form (|5H|l for I(z), i.e. A 2 = — 2slnz. Collecting both terms we finally 
have 

t2 -v 2 

A a = , ln 2 z-2slnz + 0(lnz) (A23) 
a 4 

where additional lnz terms involve A 2 (w) and t as in Eqs. (|47l I48fl . 

We proceed to identify A\(u), which determines the BG domain size, and to examine the condition A±{u c ) > zd /c- 
necessary for the appearance of the In 2 z term in l|A23[) . Eqs. (|A5I IA16|> yield for the range 2a.y/dz/c- <u<u c , 

ai{u ) = ^-^e- c ^. (A24) 
gt u c 

The definition A\{u) = ua\(u) — f a\(v)dv then leads to 

A 1 (u c ) = 2 -^- Uc e- c ^ (A25) 
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Ci is a Debye Waller factor which is small by the assumption of being well below melting, T/r <C 1. Comparing with 
(|Al5|) we identify u c » 10 4 sT 'd 2 $ / X ab 2 a 4 < 1 and Ai(u c ) « a^{g/d 3 ) 2 . 

Ai(u c ) is related to the BG domain size in the axis perpendicular to the layers Lg G = -^/c_/Ai(u c ) or in the 
ab plane i?£> G = \J c'/Ai(u c ), as identified by the q, fc cutoffs in <?i(it) (Eq. IA8fl . or by evaluating displacement 
correlations^. Hence 

f 0- 4 A Qb 3 a 3 , , , 

- — (A26) 

while Lg G = Rg G a/ (X a bV^)- These forms are valid close to decoupling [Ai(u c ) > zd /c_] or in the decoupled phase 
(z = 0). Remarkably, this result of i?BG is, up to the 10~ 4 factor, identical to that found from the dimensonal analysis 
Eq. H(j2fl with (it|,) w a 2 . We do not attempt to evaluate Rbg in the coupled phase with Ai(u c ) < zc'/c- since 
then the nonsingular phase, being neglected here, is essential for generating the proper C44. As noted above, for the 
purpose of decoupling the value of A a in the k < 1/a range for large z \zd /c- > Ai(u c )] is negligible even without 
the nonsingular phase, as seen in (|A20|) . 

The condition Ai(u c ) > zc'/c_, for the appearance of the In 2 z in Eq. (|A23(1 can be written in terms of (Eq. 
EU with ^/VjcZ = 1/AfVe, 



M! - W ^y . (A27) 



For typical BSCCO or YBCO parameters this implies a renormalized anisotropy of X^/X a b > 10 5 , i.e. fairly close 
to decoupling at z = 0. Note that Rg G /y/e can be identified as Rg G , the BG domain size in the coupled phase, as 
shown in Appendix B and section IV. 

APPENDIX B: BRAGG GLASS WITH NON-SINGULAR PHASE 

We solve here the decoupling transition with nonlinear coupling of disorder (BG effects) and with the non-singular 
phase. The E v term of Eq. (|21|l is neglected, i.e. no JG effects. This describes correctly thermal decoupling, i.e. 
the line s + 1 = 1 in Fig. 1 where JG is absent within the RSB scheme. To identify the proper 7io, we expand the 
renormalized Josephson coupling — zcos[6"(r) +0"(r)] w |z[6"(r) + 9"(r)} 2 so that with the other Gaussian terms of 
(|15|l we have 



«0 = | / ^^[Gf(q,k)\e a (q,k)\ 2 +z\0 a (q,k) + b a (q,k)\ 2 + c(q,k)q 2 \b a (q,k)f 



~Vabb a *{q,k)bP{qM- (Bl) 

K z 

Formally, one needs to perform a variation of (cos[&^(r) + 9"(r)]) = exp[— ^A a ] where 

A a =J2(\° a (<l,k) + b a (q,k)\ 2 } (B2) 

to obtain the z term in I|B1|) . This procedure was also used for decoupling in presence of columnar defects 33 . We 
proceed as in the pure case 1)540 by shifting to 

6 a (q,k)=e a (q,k)+ \ b a {q,k) (B3) 

Gf {q,k) + z 



Ho = \ [(Gj\q,k) + z)\9 a (q,k)\ 2 + G- 1 Jq,k)b a *(q,k)b f) (q : k)] 



which yields 

«„>, *) - v% + cm? + g2 + (1+ ;5,w ^ - i°" s m 

where the last term corresponds to the B 2 term of C44 (Eq. I55fl with A c replaced by A,?. Note that for q 3> 
(1 + Xatfk 2 ) 1 / 2 1 'Xf this reduces to (|A3|) with a 2 — > 0. A term corresponding to the last term of (|55|l . being ~ (A,?) -2 , 
is neglected. 
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We proceed to evaluate gx(u) with l|A8(l replaced here by 

»w - £ I>v + w*) + , 2 + (1+A l^ )/(A?)2 ]^ - *(«>]- ■ <■») 

For k > 1/a c(k)k 2 3> Ai(u) and g\(u) is Ai independent, as in Appendix A. For k < 1/a two regimes are identified, 
where the coefficient of the k 2 term in (IB5II becomes 



c+ = c (0) + z(Af) 2 = C _ + ^ 9 <1/A* 



a 4 r 



- (0) ^ M VT ?>W (B6) 

so that c + /c_ = 1 + 1/e ^> 1 with e defined in (|59|l . This reflects the significant dependence of C44 on interchanging 
the q — > and l/Af — > limits, as discussed in section IV. After the k integration we obtain (replacing I Al 1(1 

//an a 2 r /" 1/A " 1 f 1 , 3do , r— , s 

si A i / -*=+/ ~= ; = = ^«±/a/aT B7 

87rd 2 Jo 7i/a« \/c^ a/cV + Ai 

where a± = a? / (8ttc' y/c±) with a + for yj Ai(u)/c' > 1/Af and a_ for y/Ai(u)/cf < 1/Af. Hence Ai(u) = u 2 /4a± 
and Eqs. fA16l EHjl are valid in both a± regimes. Comparing (|A25|) with u 2 /4a± identifies u c » 2a±go/o' 3 and 
Ai(u c ) w a±((7o/^ 3 ) 2 - The BG scales -R^ G = y/ d / A\(u c ) are therefore 



^BG ~ ~J?2 R BG < 



10~ 4 A a b 3 a 3 

fi BG ~ — T^i R BG > V A ^ ( B8 ) 

Aa so 



so that = y/a+fct-RgQ = Rg G /y/e. The range Rg G < A,? < R% G allows for both length scales and serves as 
a crossover between the regimes in (|B8|) . The ratio R^ G = Rg G /\/e reflects the change in elastic constants, as in the 
dimensional argument of section IV. The result l|B8|) for R^ G agrees with (|A26ll in Appendix A. 
Renormalization of z requires the sum l|B2|l which is averaged with respect to Ho of (|B4(1 



A* = £ k) \ b a (q, k) + b a (q, k)\ 2 ) 

q,fe °/ yii k > + z 

= E U-V IaZ + ^T^TT—fGU^ k)} . (B9) 
^ k G f (q,k) + z G f L (q,k) + z 



The first term is w (T/r) In z and is neglected at T < t. The second term has a factor 

„2 



(BIO) 



G7 i (q,fc)+z a 2 + (1 + A Qb 2 fc2)/(A«) 2 
which for q < (1 + \ a b k^^/Xf strongly reduces the # integration, while for larger g, A a becomes 

2 

4*=/ — [1+ / ^j— r (Bll) 

where / indicates q > (1 + A a fc 2 fc 2 ) 1 / 2 /Af . For fc > 1/a, c(k)k 2 > Ai(u) provides a cutoff with the result ^4 2 = 
2sln(A c /z) as in Eq. I|A22|) . For k < 1/a the v integral term of (|B11|) becomes A\ as in l|A18|) except for a q cutoff 

in J . The k integration of (jA18|) produces a cutoff q > yfzjcl = l/(\fy/e) > (1 + A a fc 2 fc 2 ) 1/2 /Af , hence (jA19|) is 
valid. For Ai(v) < zc'/c_ 

1 /""V* d w 

/ -5-Ai(u) ~ const. (B12) 

\ Z In V z 
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while for A\(v) > zd /c- i|A21() is reproduced. The latter integration range exists if Ax(u c ) > zc'/c-, i.e. R BG < 
\J C-l z = X^y/t. Using (|B8|> we identify the condition for the appearance of the In 2 z term as 

R bg < onSGt 01 m2 z tcrm • ( B13 ) 

This is also the condition found in Appendix A (Eq. IA27I) . as well as the condition of section IV, as illustrated in 
Fig. 3, for the onset of the anharmonic regime. 



APPENDIX C: JOSEPHSON GLASS WITH NON-SINGULAR PHASE 



In this appendix we extend the solution of section III to include the nonsingular phase. In particular we identify 
the pinning length R p in the coupled phase and show that it coincides with (|62J) (with (u 2 -) k (q), up to a numerical 
prefactor. Since disorder is linearized, we do not expect to derive BG domain sizes. Also the integral I{z) is 
reconsidered. 

Consider then Eq. 1)12(1 with the pure part replaced by Eq. 1151) . The harmonic part can be written as 
Gfiq, k)\b a (q : fc) - b a (q, k)\ 2 + [c(q, k)q 2 5 af3 - s ^]b a (q : k)b?*(q, k) 

= Gf(q, k)\b a (q, k)\ 2 + d a B-l(q, k)d^(q, k) - Gj 2 (q, k)B^(q, k)b a *(q, k) (CI) 

cT(q,fc) = b a (q,k)~B^ a (q,k)Gj 1 (q,k)b a (q,k) 

q 2 

B a )}(<l' k ) = G f 1 {q,k)a{q,k)5 a} f3-s -i I 

K z 

a(q,k) = 1 + G f (q,k)c(q : k)q 2 . (C2) 
The resulting replicated Hamiltonian is 

i r c(g, fc) 2e s oQ 2 



where 



W P0 /T = £ £ fe) ^* (q , k) + \[^q 2 5 a , p - -^L^ (q, k)iP * iqi k) 

-|^E/ ^rcoB^r)-^ J d 2 rcosK(r)-P n (r)}. (C3) 



The effect of the nonsingular phase on our previous Hamiltonian Eq. 112(1 of section II is to replace c(q, k) — > 
c(q, k)/a(q, k) and s — > s /a 2 (q, k). From the definition in Eq. I(C2(I we find that a(q, fc) — 1 is either ~ q 2 or ~ fc 2 
and is small except when 

a(q, fc) - 1 = e fc < l/X ab , q < ka/X ab 

„2 

fc < l/Xab, q > ka/X ab (C4) 



IS^A 2 , 



This behavior is sufficient to eliminates the fc — > divergence of I(z) (leading to ~ \j ^pz) as shown below. 
We proceed to evaluate the fluctuations in it tr (q, fc) and identify the scale R p . From eq. ((("2(1 

(Mq, k)\ 2 ) = (2^d 2 )- 2 |J[(d Q (q, k)d a *(q, fc)) + Gf(q, k)B ja (q, k)B Ya (q, k)G a0 (q, fc)] . (C5) 

K z 

Here G a /s(q, fc) = (6 7 (q, k)b*,(q, fc)) is the solution from section III, and in the replica limit 

^G 7r (g,fc) -0 

77' 
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where terms involving Aq cancel. Hence 



kj {G f L (q,k)+c(q,k)q 2 ) 2 
2s q 2 /k 2 z Gj 2 (q,k) 



2 



Gj 2 (q,k) 



With some straightforward algebra, 



{G- f 1 {q,k)+c(q,k)<l 2 T $$<l 2 + z 



(C7) 



(|ur(q,fe)| 2 ) - (2^ 2 )- 2 | ? [.s g 2 G / (q,fc)a- 1 (q,fc) ( c(q, fc)q 2 + - - - 

s ^ c(q,fc) 2 



c(q, fc)a 2 (q, fc) V«(q> fc ) 



q 2 + 2 ] + ... (C8) 



where ... stands for terms which converge in (q, k) integration. Note the term Gj 1 (q, k)z/[Gj 1 (q, k) + z] which 
depends on the order of q — > and z — > limits; this limit dependence leads to the apparent discontinuity in C44 as 
discussed in section IV. For z^O and small q, i.e. Gj 1 (q, fc) -C 2 the first term in Eq. (|C8() dominates, leading to 

(|urM)|,)w .v^ 9<1/A * (C9) 

where C44 is from Eq. I|57|l and the condition Gj 1 (q, fc) <§C 2; is written in terms of (Eq. I56fl . The correlations at 
distance r parallel to the layers are then 

4<7 2 snT 2 r 

([u T (r) - u T (0)} 2 ) « - — r = $± . (CIO) 
a 4 c 44 c 6 ^ u p 

The last equality defines the pinning length R p where the fluctuations become of order £ 2 . This result for R p (up to 
a numerical prefactor) is the same as the one obtained from Eq. (|62|l with (u\) ~ £ 2 . 

In the decoupled phase with z — the second term in Eq. I|C8[1 dominates. To leading order in e the result is 
identical to Eq. I|C10(I except that C44 is replaced by its z = value Eq. J^HJ, i.e. the pinning length is reduced. 

Consider next the integral I(z). As noted below Eq. (|C3|) the nonsingular phase leads to the replacements 
c(q, fc) — > c(q,k)/a(q,k) and sq — > so/a 2 (q,k) so that Eq. (|46|l becomes 

j-,-, = I dq 2 dk 1 

J k 2 c(q,k)a(q,k) [c(q,k)/a(q,k)]q 2 + z 

In the range 1/a < k < n/d with a(q, fc) « 1 the q 2 term in c(q,k) amounts to a cutoff q 2 (defined below (|36|^ 
leading to 7 (z) 105}. In the range 1/A a & < fc < 1/a we have c(q, fc) = c(0)(l + 2ir\ ab 2 q 2 /a 2 k 2 ) and a(g, fc) w 1. The 
singularity in 2 which we wish to identify, is exhibited by q — > 0, hence c(q, fc) « c(0) = c_ leads to the first correction 

\ 

7!^) = ^ In 2 + consi. (C12) 

c_ 

In the range fc < 1/A a (, we have two terms 

" 1/Xab dk f X *bk/a dg 2 



i (2) - 2 r r 

7o fc 2 c_ J 



2b <r 



7 3 ( 2 ) = 2 / lMab ^ [ 1/a ( 16^A ab 2 fc 2 )2 dg 2 

Jo ^ 2 c" J\ ab k/a a V C 'V + z 

where c" = (327r 2 A a & 4 /a 4 )c_ is due to the finite effect of c(q, fc)/a(q, fc) when q > X a bk/a, k < 1/X a b- ln 7 2 (2) the q 4 
term replaces the q 2 cutoff as ak/X a b leading to 

h(z) = -^Ur\/^ + ^ m z + const - (C14) 
A q6 V V z ci 
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while 

^! :: A„,, ^ I'- 
ls smaller then the 1st term of (IC14II . We conclude then 



X nh -, :1 _ V z +2sln z 



I(z) = - 2 2 \l + 2s\n—- + const. (C16) 



The effect of the 1/ \fz term is significant, in terms of A,? (Eq. I56|l if 

Af V2\ ab 2 a ^ 4 
A ab 4d 2 ln 2 (a/rf) 

for BSCCO parameters; with bare anisotropy of \ c /\ a b ~ 50 one needs to be fairly close to the transition to have an 
effect from the 1/y/z term. Note that nonlinear coupling of disorder, i.e. BG formulation, is much more efficient in 
reducing the I(z) singularity, as shown in Appendices A and B. 
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